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Abstract 

We analyze in detail the constraints on SUSY-model parameters obtained from 
K — K mixing in the hypothesis of a splitted SUSY spectrum. FCNC contributions 
from gluino-squark-quark interactions are studied in the so-called mass insertion ap- 
proximation. We present boundaries on mass insertions and on SUSY mass scales. 
We improve previous results by including the NLO-QCD corrections to AS = 2 effec- 
tive Hamiltonian and the complete set of B-parameters for the evaluation of hadronic 
matrix elements. A full set of magic-numbers, that can be used for further analyses 
of these models, is also given. We find that the inclusion of NLO-QCD corrections 
and the B-parameters change the results obtained at LO and in the Vacuum Insertion 
Approximation by an amount of about 25 — 35%. 



1 Introduction 



It is well known that SuperSymmetry (SUSY) introduces many new sources of Flavor Chang- 
ing Neutral Currents (FCNC) which give strong constraints on the construction of extensions 
of the Standard Model (SM). 

A common feature of these models is that FCNC effects are induced by SUSY breaking 
parameters that mix different flavors. In the literature several ideas have been proposed in 
order to suppress unwanted FCNC effects. For instance, in models where SUSY breaking 
is induced by gauge interactions |] SUSY breaking parameters are flavor blind or they are 
dominated by the dilaton multiplet of string theory |2|]. Alternatively, flavor symmetries are 
used to provide either a sufficient degeneracy between the first-two generation of sfermions H 
or alignment between quark and squark mass matrices 0. 

Here we want to investigate the hypothesis that the average squark mass of the first-two 
generations is much higher then the rest of the spectrum of (s-)particles [§, || [7|, § . Through- 
out the paper we indicate the average mass of the heavy scalar squarks as M sq and the typical 
mass scale of gauginos and of the other light sparticles as m g . Small Yukawa couplings of 
the first-two generations of scalars to the Higgs doublets, together with masses of the rest of 
the supersymmetric spectrum close to the weak scale allow a natural electroweak symmetry 
breaking (EWSB). This scenario has very interesting phenomenological signatures |§ and 
can be easily realized in string theory ||. 

We consider gluino-squark-mediated FCNC contributions to AMk and e# in the neutral 
.fT-system. The effect of the most general squark mass matrix for this class of models is 
studied. In some cases further restrictions on the squark masses are required and other con- 
tributions can be more important. In particular chargino- squark- quark interactions should 
be also considered. We postpone a discussion with the inclusion of these effects to a subse- 
quent work. 

We work in the so-called mass insertion approximation flOfl . In this framework one 
chooses a basis for fermions and sfermions states where all the couplings of these particles 
to neutral gauginos are flavor diagonal and FC effects are shown by the non-diagonality of 
sfermion propagators. The pattern of flavor change, for the K-system, is given by the ratio 
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where (wi£-)ab are the off-diagonal elements of the d mass squared matrix that mixes flavor 



i, j for both left- and right-handed scalars (A,B =Left, Right), see e.g. [11]. The sfermion 
propagators are expanded as a series in terms of the 5's and the contribution of the first 
term of this expansion is considered. 

The supersymmetric flavor problem consists in building viable models in which FCNC 
are suppressed without requiring excessive fine tuning of the parameters. 

In models with a splitted spectrum of s-particles, in which the average mass of the lightest 
(trig) is in the electroweak or TeV region, two scenarios are possible: 
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for reasonable values of M sq , the suppression of FCNC requires small 8's values. Thus, 
by fixing M sq , one can find constraints on 5's, see refs. [11], [12|] and for a very recent 
NLO analysis ref. pj; 



for natural values of the 5's, say 0(1) or order of the Cabibbo angle, 0(0.22), one finds 
that the only way to get rid of unwanted FCNC effects, is by having the squarks of the 
first-two generations heavy enough. Thus, by fixing the <5's, one can find constraints on 
the minimal values for M sq . Large values of M sq , however, induce large values for the 
GUT masses of the third generation of squarks via Renormalization Group Equations 
(RGE). Consequently there can be fine-tuning problems for the Z-boson mass. We 



study this issue in sec. HI This point of view was adopted in refs. [16, 17 . 



In the past, several phenomenological analyses were carried out, which relied on some ap- 
proximations. For instance, the work of ref. |16| does not include QCD radiative corrections 
and makes use of Vacuum Insertion Approximation (VIA) for the evaluation of hadronic ma- 
trix elements. Leading order QCD corrections to the evolution of Wilson coefficients were, 



instead, considered in the papers of refs. [12, 17|. These authors found that QCD corrections 



are extremely important. For example in |T7[ the lower bound on the heavy squark mass is 
increased by roughly a factor three. 

In this work we discuss both the cited scenarios and improve previous analyses including 
the Next-to-Leading Order (NLO) QCD corrections to the most general AF = 2 effective 



Hamiltonian flTJJ and the lattice calculation of all the B-parameters appearing in the K — K 



mixing matrix elements that have been recently computed P0 |. We find this very interesting 
for several reasons. First of all we find that the inclusion of these effects leads to sizeable 
deviation from the previous computations. The results obtained using only LO-QCD correc- 
tions and Vacuum Insertion Approximation are corrected by about 25 — 35%. Furthermore, 
the uncertainties of the final result due to its dependence on the scale at which hadronic 
matrix elements and quark masses are evaluated is much reduced. 

Predictions for any model can be tested using the so-called magic numbers we provide. 
These numbers allow to obtain the coefficient functions at any low energy scale once the 
matching conditions are given at a higher energy scale. The magic-numbers will be useful, 
e.g., when a complete NLO analysis of SUSY contributions to AF = 2 processes (which 
should include also chargino exchange effects) will be implemented in the future. 

A complete NLO calculation should be comprehensive also of the 0(a s ) corrections to 
the Wilson coefficients at the scale of the SUSY masses running in the loops. So far, we miss 
this piece of information for gluino-squark contributions^. We can argue the smallness of 
these corrections from the smallness of a s at such scales. This uncertainty can be removed 
only by a direct computation. 

The paper is organized as follows. In sec. [2] we introduce the formalism concerning the 
operator basis, the Wilson coefficients and the Renormalization Group Equations (RGE). In 
sec. |3| constraints on the S's are derived. The problem of consistency of the squark spectrum 



1 The matching conditions for charged-Higgs and chargino contributions have been recently computed 
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for given entries of the 5's is considered in sec. [|. Finally our conclusions are written in sec. || 
and all the magic-numbers are given in the appendix. 



2 Effective Hamiltonian and hadronic matrix elements 

In this section we describe the framework in which the basic calculations have been per- 



formed. We follow the discussion of ref. [12| in the case M sq ^> nig. Throughout the paper 
(unless otherwise explicitly specified), we assume that the average mass of gluinos and of 
the squarks of the third generation are of the same order of magnitude. 

The three steps needed to use the Operator Product Expansion (OPE) (matching of the 
effective theory, perturbative evolution of the coefficients and evaluation of hadronic matrix 
elements) are treated in detail in the following subsections. 



2.1 Operator basis and matching of the effective theory 

In order to apply the OPE one has to calculate coefficients and operators of the effective 
theory. One first integrates out the heavy scalars of the first-two generations at the scale 
M sq . This step produces AS = 1 (of the form dggs) as well as AS = 2 operators, at the 
same order 1/MjL. When also gluinos are integrated out at rrig, AS = 1 operators generate 
AS = 2 contributions that are proportional to m~/M* g , and so can be neglected. 

The final basis of operators is: 

Q 1 = ^(1-75)5*^(1 -75) A 

Q 2 = ^(1-75)^^(1-75)^, 

Q 3 = ^(1-75)^^(1-75)^, 

Q 4 = d*{l- 7s )s a (F(1 + 7 5 )A 

Q 5 = ^(1-75)^^(1 + 75)8°, (2) 

together with operators Qi,2,3 which can be obtained from Qi,2,3 by the exchange (1 —75) <-» 
(l + 7 5 ). 

The Wilson coefficients at the matching scale M sq are (see e.g. |ll|, |12|): 
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-132f 6 (x)(5 d l2 ) LR (5' 



12) RL 



(24;c/ 6 (x) + 120/ 6 (z)) (^ 2 )ll(5 
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(3) 



where x = {m g /M sq ) 2 and 

/ 8 (aO 



6(1 + 3x) lnx + x 3 - 9x 2 - 9x + 17 
6(x- l) 5 ' 
6x(l + x) In x — x 3 — 9x 2 + 9x + 1 
3(x- l) 5 ' 



(4) 



The coefficients for the operators Qi,2,3 are the same as those of Qi,2,3 with the replacement 
L <-> i?. The authors of refs. [112], [16], [T7| use the matching coefficients directly in the limit 



x — > 0. However, we have contemplated also the extreme case of ~ M sq /2, so that we 
keep the whole expression. Of course, the value of the coefficients is the same as that of 
refs. [fL2|, 16, 17 in cases where x <C 1. 

As we said, NLO-corrections to these coefficients have not been computed yet. We assume 
they are negligible, in view of the smallness of a s (M sq ) and of the fact that similar corrections 
turned out to be rather small in the SM, the two Higgs doublet model, and for the chargino 
contribution in the constrained MSSM j22|. Our effective Hamiltonian is so affected by a 
residual renormalization scheme dependence because of the missing piece of 0(a s (M sq )) in 
the matching. 



2.2 Evolution of Wilson coefficients and running of a s 

In order to evolve the Wilson coefficients between M sq and the scale at which hadronic 
matrix elements are evaluated (yU=2 GeV), one has to account for the presence of all particles 
whose mass is intermediate between the two scales, both in the /3-function of a s and in the 
Anomalous Dimension Matrix (ADM) of the operators. 



For what concerns the former one has 25 : 



(3{a s ) = -foot, - fro?. + <D{a% (5) 
A) = -(llN c -2n f -2N c n s --n g ), (6) 
1 / 9 13N? - 3 9 4AT? - 3 3N? - 1 \ 

where N c = 3 for color SU(3) and rif is the number of fermion flavors. The terms proportional 
to rig and n q represent, respectively, the gluino and light scalar contributions. n g =l and n q =4 
when one evolves between M sq and m g and to n g = n q =0 evolving from m g to a lower mass 
scale. 



4 



In ref. JT8[ the Anomalous Dimension Matrix of the operators was computed at NLO. In 
that reference, since all SUSY particles are taken to be heavy, only loops with fermions and 
gluons were considered. This result must be modified taking into account that, from M sq to 
also the squarks of the third generation and gluinos can run in the loops. As a matter 



m 



of fact, for the K- system, light third generation squarks and gluinos can enter two-loops 
ADM only via the renormalization of the gluon propagator. An explicit calculation shows 
that the required modification consists in considering the ADM of ref. |18| as a function of 
rif + N c rig + fig/4 when one evolves between the heavy squark and gluino mass scales and 
as a function of rif below the latter scale. This substitution is no more true if the squarks 
of the first-two generations are light too. 

The value of the Wilson coefficients at the hadronic scale \i = 2 GeV where matrix 
elements are computed can then be easily calculated. Following ref. [[18] one evolves between 
two scales according to: 



C(/i) 
U[fi,M] 



N[fx] U[ji,M) N~ l [M] C(M), 



1 + 



4tt 



a s (M) 



Am), 

y (0)T /(2/3 ) 



(8) 



where ■y^ is the LO-ADM and C{jj) are the Wilson coefficients arranged in a column vector. 
This formula is correct up to the NLO. U[/j,,M] gives the LO evolution already computed 
in ref. \Y1\ while J gives the NLO corrections calculated in ref. [ IE] . J depends both on the 
number of active particles at the scale /i, and on the renormalization scheme used for its 
computation. We have used J in the same scheme used for the lattice calculation of hadronic 
matrix elements, that is the so-called Landau-RI scheme. In this way the renormalization 
scheme dependence of the final result, at the scale at which hadronic matrix elements are 
evaluated, cancels out at this perturbative order. As already been stressed, for a complete 
scheme independence of our result one should include also the NLO corrections of the Wilson 
coefficients at the high matching scale. 



We provide here the full set of Wilson coefficients at /i=2 GeV as functions of M sq 
m g (the so-called magic-numbers). We find 



and 



c*(aO = £ 

r,j=l L 
5 

CiijUg) = 

r,j=l L 



(r) a a (m- g ) m 



«" r K) Cj(m- g ), 



jr) a s (m~g) j r) a s (M sq ) (r) 



47T 



47T 




C j( M sq)- (9) 



The complete expression of a r , a' r , b\j , . . ., is given in the appendix. Eq. (Q) is useful for 
testing predictions for any model, once the two scales are fixed. The magic numbers for 
the evolution of C*i_3 are the same as the ones for the evolution of Ci_3. Eq. @ and the 
formulae of the appendix, can be used with B-parameters evaluated at ix = 2 GeV (see 
eq. (|T2"D), in order to determine the contribution to AMk and ex at NLO in QCD for any 
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model of new physics in which the new contributions with respect to the SM originate from 
the extra heavy particles. It is sufficient to compute the values of the coefficients at the 
matching scales M sq and m- g and put them in eq. (^|). 



2.3 Hadronic Matrix Elements 

The hadronic matrix elements of the operators of eq. (Q) in the Vacuum Insertion Approxi- 
mation (VIA) are: 



{K«\Q X \K*) VIA 



(K°\Q, 
(K°\Q S 



(K°\Q 4 



(K°\Q t 



K°)viA 

K°)via 
K°)via 



24 \m s + m d 



M " ' M K f K , 



M, 



K 



24 \m s + m d 



M K f K , 



ill 



K 



1 1 

24 4 Vm s + m d 



M K f K , 



K°) 



VIA 



M, 



K 



1 1 

8 12 \m s + m d 



M K f, 



K ' 



(10) 



where is the mass of the K meson and m s , m d are the masses of the s and d quarks 
respectively. An analogous definition holds for Qi t 2,3- 

Hadronic matrix elements can be evaluated non-perturbatively introducing B-parame- 
ters, defined as follows: 



(K°\QM 

{K°\QM 

(K°\QM 
(K°\QM 



K°) = -M K fjkBi(ji), 



K°) 
K°) 
K°) 
K°) 



M, 



K 



24 \m s (fi) + m d (fi) j 



M 



K 



24 lm s (ju) + m d (ji) 



M K f*BM, 



M K f 2 K B 3 (fi), 



M 



K 



4 \m 8 (n) + m d (n) 



M, 



K 



12 \m s (fi) + m d (fi) 



M K f 2 K B^), 



M K f 2 K B 5 (fi), 



where Qi(fi) are the operators renormalized at the scale /i. The B-parameters for Qi^^in) 
are the same as those of <3i,2,3(/-0- 

In the computation of Bi for the operators 2-5, smaller contributions of higher order in 
chiral expansion, coming from axial current, have been neglected. A detailed explanation of 
the reasons of this approximation can be found in ref. |2(J. The definition of B-parameters 
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Constants 


Values 


a em (M z ) 


1/127.88 


a s (M z ) 


0.119 


M K 


497.67 MeV 


fx 


159.8 MeV 


m d (2 GeV) 


7 MeV 


m s (2 GeV) 


125 MeV 


m c 


1.3 GeV 


m b 


4.3 GeV 


m t 


175 GeV 


sin 2 6 W (M Z ) 


0.23124 



Table 1: Constants used for phenomenological analysis. 



in eq. ((TTJ) takes explicitly into account this approximation and using it the low scale (/i) 
dependence of the final result is explicitly canceled in the product of coefficient functions 
and hadronic matrix elements. 

The B-parameter of the first operator is usually addressed as B^ and has been extensively 
studied on the lattice and used in many phenomenological applications (see, e.g. |23], ^|). 
We have considered its world average [23||. The other Bi have been taken from ref. p(J (for 



another determination of these Bi, calculated with perturbative renormalization see ref. 

All the B-parameters are evaluated at a scale of 2 GeV in the LRI renormalization scheme: 



B 1 ( fi = 


2 GeV) 


= 0.60 ±0.06, 


B 2 (» = 


2 GeV) 


= 0.66 ±0.04, 


B 3 (fi = 


2 GeV) 


= 1.05 ±0.12, 


B^ = 


2 GeV) 


= 1.03 ±0.06, 


B S (ji = 


2 GeV) 


= 0.73 ±0.10. 



(12) 

So far in the literature all phenomenological analyses on this subject have used the VIA 
and have computed Wilson coefficients and quark masses at a scale variable between 0.5-1 
GeV. We will see this represents in some cases quite a rough approximation. 

Finally we give in table ([!]) all the numerical values of the physical constants we have 
considered. All coupling constants and sin 2 6 W (M Z ) are meant in the MS'-scheme [^6 . 
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3 Constraints on the £'s 



We are ready to provide a set of constraints on SUSY variables coming from the Kl — K$ 
mass difference, AMk and the CP violating parameter ex defined as 



AM K = 2Re(K°\H eS \K 



0\ 



£ *' = v^ lm{KW ° } - (13) 

The parameter space is composed of two real and four complex entries, that is M sq , rrig and 

{$12) LL, (812) LR, ($12)RL, (812) RR- 

Neglecting interference among different SUSY contributions, we give upper bounds on 
the 5's, at fixed values of M sq and rrig, with the condition M sq > rrig. In this way one gets 
a set of constraints on individual S's. Indeed, since we are interested in model independent 
constraints, it is meaningful to study the interference of cancellation effects only in specific 
models. 

The physical condition used to get the bounds on the <5's is that the SUSY contribution 
(proportional to each single 5) plus the SM contribution to AMk and ex do not exceed the 
experimental value of these quantities. For what concerns the SM contribution to AM^, 
we assume that the values of the CKM elements V c d and V cs are unaffected by SUSY. This 
implies the (very reasonable) hypothesis that SUSY does not correct significantly tree level 
weak decays. The value of the SM contribution to ex, instead, depends on the phase of the 
CKM matrix. This phase can be largely affected by unknown SUSY corrections and can be 
treated as a free parameter. We put the CKM phase to zero so that the experimental value 
of ex is completely determined by SUSY. Finally, to be even more conservative, we subtract 
one standard deviation to the values of the B-parameters. 

The final results are shown in tabs. |2|-|7] for gluino masses of 250, 500, 1000 GeV. We 
consider the heavy squark masses expected in some common models (see e.g. H || [7]. || ) . 

The constraints that come from the four possible insertions of the 5's are presented: in 
the first and second rows only terms proportional respectively to (6( 2 )ll and (Sf 2 )LR are 
considered; in the last two rows the contribution of operators with opposite chirality, RR 
and RL, is also evaluated by assuming (Sf 2 )LR = (812) rl and (Sf 2 )LL = ($12) rr- 

In each column of the table we show the bounds on the 5's in the various approximations 
that one can use for their determination: without QCD correction and in VIA, with LO- 
QCD corrections and in VIA, with LO-QCD corrections and with lattice B-parameters and, 
eventually, with NLO-QCD corrections and lattice B-parameters. Comparing the values 
of our constraints at LO-VIA with those found from the authors of ref. |12| we find some 
differences. The reason is twofold. On the one hand they do not consider the SM contribution 
to AMk and on the other they evaluate the hadronic matrix elements at a scale jl such that 
a s (p) = 1. This latter choice may be questionable because at this scale strong interactions 
break perturbation theory. 

The combination of B-parameters and NLO-QCD corrections change the LO-VIA results 
by about 25 — 35%. As expected flI2|, the tightest constraints are for the cases (Sf 2 )LL — 
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($12)rr an d ($12) lr = (8f 2 )RL- In these cases the coefficients proportional to (Sf 2 ) rr, 
(^12) £#(^12) -fix dominate the others. 

We have checked that the uncertainties of the results due to higher perturbative orders, 
are sizeable, being, in some cases up to 10%. 



4 Constraints on squarks spectrum 

In this section, following the discussion of ref . [ IE| , we provide a different kind of constraints. 

For fixed values of the 5's and of the average light sparticle mass, rrig, it is possible to 
calculate the minimum value of M sq necessary to suppress the FCNC at an experimentally 
acceptable level. Here we give constraints on M sq and we discuss about their consistency. 
Using Renormalization Group Equations, one finds that a too large M sq can drive to zero or 
negative values the average mass of the third generation of sfermions, mj, at the TeV scale 
(m^(~ ITeV)). To circumvent this problem, a minimum value for mj(/^GUT) at the GUT 
scale has to be chosen. If m^Gux) is too high (say more then 3-4 TeV), however, a too large 
fine-tuning of the SUSY parameters is required in order to account for the observed mass 
of the Z-boson and severe naturalness problems arise |TJ], [TJ| . This problem was studied in 



refs. [ 1~6[ , |17 |. 



One obtains constraints about the consistency of models with a splitted mass spectrum 
following three steps: 

• determining the minimum value of M sq necessary to suppress FCNC. This is discussed 



in subsec. PO; 



computing the maximum value of M sq allowed by positiveness of light scalar masses 
and fine-tuning. More about this in subsec. 



combining the previous two results one can determine regions of allowable values of 
M sq that satisfy both the requests of the previous points. We comment about that in 
subsec. Ol. 



4.1 Minimum values for heavy squark mass 

In order to obtain constraints on M sq one has to specify a value for the 5's. We consider the 

cases 
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Figure 1: Lower bounds on M sq from AM^ with various approximations for the case I with 
K, = 0.22. In this case the larger correction to LO-VIA come from the B-parameters. 



where /C can take the values (1, 0.22, 0.05). We have chosen these entries to leave aside 
possible accidental cancellations. The cases in which K.—1, are, of course, extreme cases: 
one may wonder about the consistency of mass insertion approximation as the neglected 
terms are of order 0(1). However these cases have already been studied in the literature, 
|l7|), and we report them for completeness. The results so obtained just give 



16 



(see e.g. 

an estimate of the mass scales that are involved and can be trusted if other corrections do 
not provide accidental cancellations. This can be checked only by a direct calculation. 

The assumptions made for the SM contribution and the B-parameters are the same as 
in section |3|. 

In order to monitor the effect of the different corrections on the final result we show in 
figs. [1], H the lower bound obtained for the cases I and III with /C = 0.22 (the other cases 
give similar results). As we see, B-parameters and NLO-QCD corrections play a significant 
role in the final computation and the correction they provide with respect to the LO-VIA 
results are of the order of (25 — 35)%. In particular, in case I, fig. [I], B-parameters provide 
the most important corrections with respect to LO-VIA results. In case III, fig. |2|, instead, 
corrections to LO-VIA results are dominated by the NLO-QCD perturbative contributions. 

The case (Sf 2 )LL = (^12) rr was al so considered in refs. [T7|] . The differences, at LO 



and without B-parameters, among our result and the ones of refs. [jig , |17f come from our 
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Figure 2: Lower bounds on M sq from AMk with various approximations for the case III 
with K, = 0.22. In this case the larger correction to LO-VIA come from NLO perturbative 
corrections. 



inclusion of the SM contribution, from the value of the strange quark mass and from the 
scale at which hadronic matrix elements are evaluated. We agree with them for the same 
choice of parameters. 

Notice that, if the imaginary parts of the <5's are of the same order of their real parts, 
there are much stronger constraints coming from ex than from AMk (namely by a factor 
~ 7.7). To be conservative, we consider in this section only constraints coming from the real 
parts of the <5's. 

The final results are shown by the (colored) continuous lines in figs. [], f|, ^. The minimum 
value of M sq depends strongly both on K and on the case one considers (I, II, III or IV, 
see eq. |TJ). Notice that (Sf 2 ) LR , (Si 2 )rl, (which enter the cases II and IV), are "naturally" 



small in the MSSM. However, since we would like to do a model-independent analysis, we 
have made no particular assumption on them. In all graphs the strongest constraints come 
from the case (Sf 2 )iR = (Sf 2 ) rl ^ 0. Much lower constraints are generally obtained in cases 
I and II. In fig. || the case I has not been drawn since no constraint can be derived. 
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Figure 3: The full (colored) lines give the lower bounds on M sq necessary to suppress FCNC 
and with JC=1 for the various cases. An upper bound on M sq is derived in order to satisfy 
fine-tuning requirements and it is shown by the dashed line. The two kind on constraints are 
not compatible in this case. 



4.2 RGE for the masses of the third generation of scalars 

It is well known that large values of M sq can drive the mass of the third generation of scalars 
to negative values, via RGE ||16|| . Let us consider the two- loop RGE's for the mass, mj of the 



third generation of scalars, /. In the DR scheme (see e.g. ref. [27]), with two generations 
of heavy scalars, one has 

= -^Il^cRmlU^ + j^Yl^clM^ (15) 

where C{ is the Casimir factor for / in the SU(5) normalization, the sums are over the 
gauge groups SU(3), SU(2), U(l) and mc denotes the gaugino masses. In eq. (^), Yukawa 
couplings are neglected: these couplings drive the light masses to even lower values and 
so, in this respect, our choice is a conservative one. Moreover, the introduction of Yukawa 
interactions requires further assumptions on SUSY parameters (see e.g. [|I7J]) that we do not 
discuss in this paper. 

The solution of eq. (|TJ) between a Grand Unification (GUT) scale /igut ~ 2 • 10 16 GeV, 
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0.50 



1.00 1.50 
rrig [TeV] 



2.00 



Figure 4: The same as fig. [| with JC=0.22. Cases I and II are now compatible with fine-tuning 
requirements. 



and /i ~ 1 TeV can be easily written as 



w/Ogut) - E 



16 

4^ 



^(//gut) - a i (M sg )lc/ M 2 sq + 



i Pi 



E4[(™g) 2 (m 

i Pi 



(m^)i(^) 



(16) 



where are the /3-functions LO coefficients of the z-th gauge coupling. In eq. ([IBD we have 
considered a common gaugino mass, mc at the GUT scale, while for what concerns the 
couplings we have evolved them starting backward from \x = Mz- Note that in eq. (16) the 
contribution of the heavy scalar s has been decoupled at M sq . 

Eq. (|16"D can be used in order to derive consistency constraints on the values of M sq and 
to|-(/4gut) once the values of m\(fi) and of m^-(/i GUT ) are fixed. The latter can be determined 
according to the following requirements. First, m 2 j(fi) must be at least positive, such as to 
leave color and electric symmetries unbroken. The value of ^^(/igut) determines the amount 
of fine-tuning necessary in order to achieve the electroweak symmetry breaking. Following 
ref. |14j the necessary fine-tuning scales approximately as 10% x (0.3 TeV/mQ 3 (/xcuT)) 2 for 



the squark doublet of the third generation Q 3 . We have calculated the constraints on M. 
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0.50 1.00 1.50 2.00 

rrig [TeV] 

Figure 5: The same as fig. [| with JC=0.05. Case I is not drawn since no lower bound on 
M sq can be obtained in this case. Cases II and III are now compatible with fine-tuning 
requirements. 
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and m^(/iGUT) coming from eq. ([Tj]) in the case f = Q3 choosing for ^^ 3 (/icuT) the value 
of (3.5 TeV) 2 . The latter choice corresponds to a fine-tuning of more than 0.1%. 

At fixed values of m 2 ^ (//) and (/xgut) (which depend on M sq and m G ) one can plot 
the upper value of M sq as function of mg. The result is the (black) dashed line of figs. [| £|. ||. 
One finds that M sq can not be much larger than about 25 TeV. Of course this is just an 
estimate of this limiting value. The inclusion of Yukawa couplings, of more severe fine-tuning 
requirements and of other effects can only lower this limit. 



4.3 Final remarks 

In figs. H|, |5] we combine the constraints derived in the two previous subsections. These 
figures (together with tables tabs. 0-(7|) suggest that also models with a splitted mass spectrum 
need further assumptions to be phenomenologically viable, e.g. one has to introduce flavor 
symmetry or dynamical generation of degenerate scalar masses ITS . 



In particular, without these further hypotheses, most of the cases which we have consid- 
ered, face fine tuning problems. In particular, values of /C ~ are hardly acceptable. 
Although /C ~ (9(0.22) and /C ~ 0(0.05) have better chances they must be treated carefully. 



5 Conclusions 

In this work we analyze in detail the constraints on SUSY-models parameters coming from 
K — K oscillations in the hypothesis of a splitted SUSY spectrum. FCNC contributions 
coming from gluino-squark- quark interactions, working in the so-called mass insertion ap- 
proximation, have been considered. We provide boundaries on mass insertions and on SUSY 
mass scales and we discus their consistency. Previous results including NLO-QCD correc- 
tions to AS = 2 effective Hamiltonian, B-parameters for the evaluation of hadronic matrix 
elements have been improved. A full set of magic-numbers is provided, that can be used for 
further analyses. 

We have discussed the residual uncertainty of our results coming from our ignorance of 
NLO-QCD corrections to the matching coefficients. 

Our analysis confirms that a splitted sparticle mass spectrum does not explain easily 
FCNC suppression without some amount of fine-tuning. These problems can be solved only 
if further assumptions in these kind of models are made, e.g. flavor symmetry or dynamical 
generation of degenerate scalar masses ||16|| . 

In order to perform a complete analysis of SUSY-FCNC effects chargino contributions 
should be included. It is also interesting to extend this kind of analysis to AB = 2 pro- 
cesses, once the calculation of B-parameters for the B — B system parameters (which is in 
progress [0) will be completed. 
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Appendix 



We give here the numerical values for the magic numbers of eq. 
entries are shown: 



Only the non- vanishing 



h (r) 
°11 
h (r) 
°22 
h (r) 
°23 
h (r) 
°32 
h (r) 
°33 
h (r) 

b (r) 
"45 
h (r) 
"54 
Jr) 
°55 



Ar) 
Jll 
Ar) 
J 22 
Ar) 
J 23 
Ar) 
J 32 
Ar) 
J 33 
Ar) 
J 44 
Ar) 
J 45 
Ar) 
J 54 

Ar) 
J 55 



0.29,-1.1,0.14,-0.69,0.79) 
0.46,-1.8,0.23,-1.1,1.3) 



5,0,0,0,0) 

0.0048,1.1,0,0) 

-0.0073,0,0,0) 

-0.23,0.47,0,0) 

0.34,0,0,0) 

0,0,0.52,-0.017) 

0,0,0.99,-2.2) 

0,0,-0.00051,0.020) 

0,0,-0.00096,2.5) 



(r) 

(r) 
c 22 

Jr) 
c 23 

(r) 
c 32 
>) 
c 33 

(r) 
C 44 

(r) 

c 45 

>) 
c 54 

°55 



Jr) _ 
"11 _ 


(1 


0,0,0,0,0) 


Jr) - 
e ll — 


(-1.5,0,0,0,0) 


d {T) ~ 
"22 — 


(0 


0,1.0,0,0) 


(r) 
e 22 — 


(0,0.70,-4.9,0,0) 


"23 — 


(0 


0,0,0,0) 


Jr) _ 
e 23 — 


(0,-1.1,0,0,0) 


d ir) ~ 
"32 — 


(0 


-0.67,0.67,0,0) 


(r) 
e 32 — 


(0,-27,-9.0,0,0) 


"33 _ 


(0 


1.0,0,0,0) 


Jr) _ 
e 33 _ 


(0,40,0,0,0) 


Ar) _ 
U44 — 


(0 


0,0,1.0,-0.015) 


(r) 


(0,0,0,22,0.46) 


"45 — 


(0 


0,0,1.9,-1.9) 


Jr) _ 
e 45 — 


(0,0,0,41,58) 


"54 — 


(0 


0,0,-0.0081,0.0081) 


(r) 
e 54 — 


(0,0,0,0.15,-0.12) 


"55 — 


(0 


0,0,-0.015,1.0) 


(r) 
e 55 _ 


(0,0,0,0.29,-15) 



5,0,0,0,0) 

0,4.2,0,0) 

0,1.1,0,0) 

33,2.8,0,0) 

-41,0.70,0,0) 

0,0,-23,0.40) 

0,0,-72,-28) 

0,0,0.18,-0.21) 

0,0,0.57,15) 



-3.4,0,0,0,0) 

0,-0.12,2.8,0,0) 

0,0.12,1.2,0,0) 

0,6.3,2.0,0,0) 

0,-6.2,0.88,0,0) 

0,0,0,-5.7,1.8) 

0,0,0,-26,-27) 

0,0,0,0.0086,-0.77) 

0,0,0,0.040,12) 
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M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 




M^ 2 )Il\ 


2 
5 

10 


3.1 x 10~ 2 
7.5 x 10~ 2 
0.15 


3.6 x 10~ 2 
8.8 x 10~ 2 
0.18 


4.9 x 10~ 2 
0.12 
0.25 


4.9 x 10~ 2 
0.12 
0.24 




MSf 2 ) 2 LR \ 


2 

5 

10 


2.1 x 10~ 2 
9.8 x 10~ 2 
0.34 


1.4 x 10" 2 

6.5 x 10~ 2 
0.22 


1.8 x 10~ 2 
8.2 x 10" 2 
0.28 


1.6 x 10~ 2 
7.2 x 10~ 2 
0.25 


^/\Re(5f 2 ) LR = Re(5f 2 ) RL \ 


2 
5 

10 


6.6 x 10~ 3 
1.5 x 10" 2 
3.0 x 10~ 2 


3.5 x 10~ 3 
7.7 x 10~ 3 
1.5 x 10- 2 


3.5 x 10~ 3 
7.9 x 10~ 3 
1.5 x 10~ 2 


2.8 x 10~ 3 
6.4 x 10~ 3 
1.2 x 10~ 2 


y/lMS&li. = M^f 2 ) RR \ 


2 
5 

10 


1.1 x 10" 2 
4.1 x 10- 2 
0.10 


5.2 x 10- 3 
1.6 x 10~ 2 
3.6 x 10~ 2 


5.1 x 10" 3 
1.6 x 10~ 2 
3.4 x 10~ 2 


4.1 x 10~ 3 
1.3 x 10~ 2 
2.7 x 10~ 2 



Table 2: Limits on Re(5f 2 )AB from AM K with gaugino masses of 250 GeV. 



M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 


A 


Re(^)LI 


2 
5 

10 


3.8 x 10~ 2 
8.1 x 10~ 2 
0.16 


4.5 x 10~ 2 

9.6 x 10~ 2 
0.19 


6.1 x 10~ 2 
0.13 
0.26 


6.1 x 10~ 2 
0.13 
0.26 


J 


Msf 2 )l R \ 


2 

5 

10 


1.6 x 10~ 2 
6.3 x 10~ 2 
0.21 


1.1 x 10~ 2 

4.2 x 10~ 2 
0.14 


1.3 x 10~ 2 
5.3 x 10~ 2 
0.17 


1.2 x 10~ 2 
4.7 x 10~ 2 
0.15 


^\Re(5 d 2 ) LR = Re(5 d 2 ) RL \ 


2 

5 

10 


9.6 x 10 -3 

1.7 x 10~ 2 
3.2 x 10~ 2 


4.6 x 10~ 3 

8.5 x 10~ 3 

1.6 x 10~ 2 


4.5 x 10~ 3 
8.7 x 10~ 3 

1.6 x 10~ 2 


3.6 x 10~ 3 
7.0 x 10~ 3 
1.3 x 10~ 2 


\/\Ms d u)lL = Msf 2 ) 2 RR \ 


2 
5 

10 


8.6 x 10~ 3 
3.2 x 10~ 2 
8.8 x 10~ 2 


4.4 x 10~ 3 
1.4 x 10~ 2 
3.4 x 10~ 2 


4.4 x 10~ 3 
1.4 x 10~ 2 
3.2 x 10~ 2 


3.6 x 10~ 3 
1.1 x 10~ 2 
2.6 x 10~ 2 



Table 3: Limits on Re(5f 2 )AB from AM K with gaugino masses of 500 GeV. 
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M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 




Re(5f 2 )i L | 


2 
5 

10 


5.9 x 10~ 2 
9.6 x 10~ 2 
0.17 


6.9 x 10~ 2 
0.11 
0.21 


9.4 x 10~ 2 
0.15 
0.28 


9.3 x 10~ 2 
0.15 
0.28 




MSf 2 ) 2 LR \ 


2 

5 

10 


1.4 x 10~ 2 
4.6 x 10~ 2 
0.14 


9.7 x 10~ 3 
3.0 x 10~ 2 

8.8 x 10~ 2 


1.2 x 10~ 2 
3.9 x 10" 2 
0.11 


1.1 x 10~ 2 
3.4 x 10~ 2 
9.8 x 10~ 2 


^/\Re(5f 2 ) LR = Re(5f 2 ) RL \ 


2 
5 

10 


4.2 x 10~ 2 
2.2 x 10~ 2 
3.6 x 10~ 2 


9.8 x 10~ 3 
1.1 x 10~ 2 
1.8 x 10- 2 


7.8 x 10~ 3 
1.1 x 10~ 2 
1.8 x 10~ 2 


6.0 x 10~ 3 
8.5 x 10~ 3 
1.4 x 10~ 2 


y/lMS&li. = M^) RR \ 


2 
5 

10 


8.0 x 10~ 3 
2.5 x 10~ 2 
6.8 x 10~ 2 


4.2 x 10- 3 
1.2 x 10~ 2 
2.9 x 10~ 2 


4.3 x 10- 3 
1.2 x 10~ 2 
2.9 x 10~ 2 


3.5 x 10~ 3 
9.7 x 10~ 3 
2.3 x 10~ 3 



Table 4: Limits on Re(5f 2 )AB from AM K with gaugino masses of 1000 GeV. 



M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 


A 


M*? 2 )Ll 


2 
5 

10 


4.0 x 10 -3 
9.7 x 10 -3 
2.0 x 10~ 2 


4.7 x 10~ 3 
1.2 x 10~ 2 
2.4 x 10" 2 


6.4 x 10~ 3 
1.6 x 10~ 2 
3.2 x 10~ 2 


6.4 x 10~ 3 
1.6 x 10" 2 
3.2 x 10" 2 




Im« 2 )Ll 


2 

5 

10 


2.7 x 10 -3 
1.3 x 10~ 2 
4.5 x 10" 2 


1.8 x 10~ 3 
8.4 x 10~ 3 

2.9 x 10~ 2 


2.3 x 10~ 3 
1.1 x 10~ 2 
3.7 x 10~ 2 


2.1 x 10~ 3 
9.4 x 10~ 3 

3.2 x 10~ 2 


J\Im(5f 2 ) LR = Im(Sf 2 ) RL \ 


2 

5 

10 


8.6 x 10" 4 
2.0 x 10~ 3 
3.9 x 10~ 3 


4.5 x 10" 4 
1.0 x 10" 3 
2.0 x 10~ 3 


4.6 x 10" 4 
1.0 x 10" 3 
2.0 x 10~ 3 


3.7 x 10" 4 
8.3 x 10~ 4 
1.6 x 10~ 3 


v /|Im( ( 5f 2 ) 2 L = lm(5 d 12 ) RR \ 


2 
5 

10 


1.4 x 10~ 3 
5.4 x 10- 3 
1.4 x 10~ 2 


6.7 x 10~ 4 
2.1 x 10" 3 
4.7 x 10" 3 


6.6 x 10~ 4 
2.0 x 10" 3 
4.5 x 10" 3 


5.4 x 10~ 4 
1.6 x 10~ 3 
3.6 x 10~ 3 



Table 5: Limits on lm(8f 2 )AB from e K with gaugino masses of 250 GeV. 
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M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 




M*? 2 )LI 


2 
5 

10 


5.0 x 10~ 3 

1.1 x 10" 2 
2.1 x 10~ 2 


5.9 x 10~ 3 
1.3 x 10~ 2 
2.5 x 10~ 2 


8.0 x 10~ 3 
1.7 x 10~ 2 
3.4 x 10~ 2 


7.9 x 10~ 3 
1.7 x 10~ 2 
3.3 x 10~ 2 


A 


Im(<5f 2 )L| 


2 

5 

10 


2.0 x 10~ 3 
8.3 x 10~ 3 
2.7 x 10~ 2 


1.4 x 10" 3 

5.5 x 10~ 3 
1.8 x 10~ 2 


1.8 x 10~ 3 
7.0 x 10~ 3 
2.2 x 10~ 2 


1.6 x 10~ 3 
6.2 x 10~ 3 
2.0 x 10~ 2 


y/\Tm(Sf 2 ) LR = lm(5f 2 ) RL \ 


2 
5 

10 


1.3 x 10~ 3 
2.2 x 10 -3 
4.2 x 10~ 3 


6.0 x 10~ 4 

1.1 x 10~ 3 
2.1 x 10- 3 


5.9 x 10~ 4 
1.1 x 10~ 3 
2.1 x 10" 3 


4.7 x 10" 4 
9.1 x 10~ 4 
1.7 x 10~ 3 


y/\M^) 2 LL = MSf2) RR \ 


2 
5 

10 


1.1 x 10~ 3 

4.2 x 10~ 3 
1.1 x 10~ 2 


5.8 x 10" 4 

1.9 x 10~ 3 
4.4 x 10~ 3 


5.8 x 10" 4 
1.8 x 10~ 3 
4.2 x 10~ 3 


4.7 x 10~ 4 
1.5 x 10~ 3 
3.4 x 10~ 3 



Table 6: Limits on lm(8f 2 )AB from t K with gaugino masses of 500 GeV. 



M sq [TeV] No-QCD, VIA LO-VIA 


LO, Bi 


NLO, Bi 


A 


Im(5f 2 )L| 


2 
5 

10 


7.7 x 10~ 3 
1.3 x 10~ 2 
2.3 x 10~ 2 


9.0 x 10~ 3 
1.5 x 10~ 2 
2.7 x 10~ 2 


1.2 x 10~ 2 
2.0 x 10~ 2 
3.7 x 10~ 2 


1.2 x 10~ 2 
2.0 x 10" 2 
3.6 x 10" 2 




Im« 2 )Ll 


2 

5 

10 


1.8 x 10 -3 
6.0 x 10~ 3 
1.8 x 10" 2 


1.3 x 10~ 3 

4.0 x 10~ 3 

1.1 x 10~ 2 


1.6 x 10~ 3 
5.0 x 10~ 3 
1.5 x 10~ 2 


1.4 x 10~ 3 

4.5 x 10~ 3 
1.3 x 10~ 2 


v /|Im(^ 2 ) LR = Im(6f 2 ) RL \ 


2 

5 

10 


5.5 x 10 -3 
2.9 x 10~ 3 
4.7 x 10~ 3 


1.3 x 10~ 3 

1.4 x 10" 3 
2.3 x 10~ 3 


1.0 x 10~ 3 
1.4 x 10" 3 
2.3 x 10~ 3 


7.8 x 10" 4 
1.1 x 10~ 3 

1.9 x 10~ 3 


v /|Im( ( 5f 2 ) 2 L = lm(5 d 12 ) RR \ 


2 
5 

10 


1.0 x 10~ 3 
3.3 x 10~ 3 
8.9 x 10 -3 


5.5 x 10~ 4 

1.6 x 10" 3 
3.8 x 10" 3 


5.6 x 10~ 4 

1.6 x 10" 3 

3.7 x 10" 3 


4.6 x 10~ 4 
1.3 x 10~ 3 
3.0 x 10~ 3 



Table 7: Limits on ha(Sf 2 )AB from 6k with gaugino masses of 1000 GeV. 
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